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Abstract 



We consider the stochastic background of gravity waves produced by first- 
order cosmological phase transitions from two types of sources: colliding bub- 
bles and hydrodynamic turbulence. First we discuss the fiuid mechanics of 
relativistic spherical combustion. We then numerically collide many bubbles 
expanding at a velocity v and calculate the resulting spectrum of gravitational 
radiation in the linearized gravity approximation. Our results are expressed 
as simple functions of the mean bubble separation, the bubble expansion ve- 
locity, the latent heat, and the efficiency of converting latent heat to kinetic 
energy of the bubble walls. A first-order phase transition is also likely to 
excite a Kolmogoroff spectrum of turbulence. We estimate the gravity waves 
produced by such a spectrum of turbulence and find that the characteristic 
amplitude of gravity waves produced is comparable to that from bubble colli- 
sions. Finally, we apply these results to the electroweak transition. Using the 
one-loop effective potential for the minimal electroweak model, the character- 
istic amplitude of gravity waves produced is /i ~ 1.5 X 10~^^ at a characteristic 
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frequency of 4.1 X 10"'^ Hz corresponding to f2 ~ 10~^^ in gravity waves, far 
too small for detection. Gravity waves from more strongly first-order phase 
transitions, including the electroweak transition in non-minimal models, have 
better prospects for detection, though probably not by LIGO. 

04.30. +x, 98.70.Vc, 98.80. Cq 
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I. INTRODUCTION 



First-order phase transitions in the early Universe can be potent sources of gravitational 
radiation [1,2]. In a recent series of papers we have calculated the radiation emitted by 
colliding vacuum bubbles and obtained useful approximations to the bubble dynamics, and 
applied these results to very strongly first-order phase transitions which occur through nu- 
cleation and percolation of vacuum bubbles [3-5]. In this paper, we extend these results to 
more weakly first-order phase transitions which occur in a thermal environment, and apply 
our results to the electroweak phase transition. 

In a first-order phase transition, the Universe starts in a metastable high-temperature 
phase (the "symmetric" phase) and converts to a stable low-temperature (the "broken- 
symmetry") phase. The transition proceeds via nucleation of bubbles of the low-temperature 
phase within the high-temperature phase; these bubbles then expand and merge, leaving the 
Universe in the broken-symmetry phase. 

Previously, we considered vacuum transitions, in which the only component of the Uni- 
verse is a scalar field. In this case true-vacuum bubbles are nucleated through quantum 
tunneling [6]. The dynamics of these bubbles is comparatively simple: once the bubbles are 
nucleated, the scalar field simply evolves according to the Klein-Gordon equation. Bubbles 
that are larger than a critical size begin to expand and rapidly approach velocities near 
the speed of light. All of the liberated vacuum energy goes into accelerating the bubble 
walls, which become progressively thinner and more energetic as the bubbles expand. These 
high velocities and large energy densities provide the necessary conditions for generating 
large amounts of gravitational radiation, and the resulting radiation spectrum depends very 
simply on the natural length and energy scales of the problem. 

For a thermal transition, the problem is more complex. Nucleation of bubbles of the 
low-temperature phase occurs through quantum tunneling and thermal fluctuations. How- 
ever, the evolution of these bubbles is not driven simply by scalar-fleld evolution. Instead, 
it depends on interactions of the bubble wall with the plasma and on the resulting fluid 
dynamics. Part of the latent heat released in the transition raises the plasma temperature, 
while another fraction of the latent heat is converted to bulk motions of the fluid. If the 
Reynolds number of the universe at the phase transition is large enough, then bubble motion 
produces turbulence in the plasma. 

In this paper, we perform detailed calculations of the gravitational radiation produced by 
the collision of spherical combustion bubbles expanding at a velocity f, using the linearized 
gravity approximation. The resulting spectra are simply expressed in terms of f, the log- 
arithmic derivative of the bubble-nucleation rate /3, the ratio of vacuum to thermal energy 
density a, and an efficiency factor k giving the fraction of vacuum energy which goes into 
kinetic energy of bulk motions of the fluid, as opposed to heating. As discussed below, the 
theory of relativistic combustion gives k and f as a function of a, which in effect, measures 
the degree of supercooling (i.e., how strongly flrst order the phase transition is). 

Combustion occurs via two distinct modes: detonation and deflagration. Roughly, deto- 
nations occur when the phase boundary propagates faster than the speed of sound, while for 
deflagrations the phase boundary propagates slower than the sound speed. We show that 
the bubble collisions in phase transitions proceeding via detonation will produce substantial 
gravitational radiation. In contrast, production of gravitational radiation from collisions of 
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deflagration bubbles should be small, because the bubble velocities are small (subsonic). It 
has recently been argued that detonation is the only stable mode of combustion for a cosmo- 
logical phase transition, and that a transition which begins via deflagration rapidly becomes 
unstable and converts to detonation [7]. For these reasons we mainly focus on gravity waves 
produced by detonation bubble collisions. 

Both modes of combustion can stir up turbulence on scales comparable to the bubble 
size. We estimate the gravity waves produced by a fully developed Kolmogoroff spectrum 
of turbulence through simple dimensional arguments, and flnd that the amplitude of the 
spectrum is comparable to that from bubble collisions. This source will generate gravity 
waves in addition to those produced by the actual bubble collisions. We note that our 
estimates are completely general, and apply to any injection of energy in the early Universe 
on a large length scale. 

Section II discusses the relevant results from relativistic combustion theory. We review 
the solutions to the hydrodynamic equations of motion for spherically symmetric detona- 
tion bubbles [8] and derive relationships between bubble-expansion velocity, bubble kinetic 
energy, latent heat, and temperature. We also discuss the solutions for spherical, relativis- 
tic deflagration bubbles, which have not been previously addressed, and briefly compare 
with the hydrodynamics associated with planar combustion [9-11]. In Section III, we re- 
view gravity-wave formalism used for our calculations. The calculation of the gravitational 
radiation produced by many colliding bubbles is made tractable through the envelope ap- 
proximation [5]; we discuss the applicability of this approximation to combustion bubbles. 
Then we numerically calculate the radiation spectra for the collision of many bubbles in 
terms of their expansion velocity and kinetic energy, which are related to parameters of 
the phase transition in Section II. Estimates of gravity waves from turbulence conclude the 
section. Section IV contains the necessary formulas to propagate the generated spectrum of 
gravity waves to the present time. As a sample application, we derive the gravitational ra- 
diation produced by the electroweak transition, using the one-loop effective potential of the 
minimal standard model. We conclude by briefly considering detection prospects, especially 
for more strongly flrst-order phase transitions. Appendix A analyzes spherical relativistic 
deflagration bubbles, and in Appendix B a model effective potential is analyzed and applied 
to the electroweak transition. 

II. FLUID FLOW IN SPHERICAL COMBUSTION 

In order to calculate the spectrum of gravitational radiation from colliding bubbles, we 
need to know the spatial components of the traceless part of the stress-energy tensor, Tij. 
For a relativistic fluid, this is simply Tij = w^f'^ViVj, where w = e -\-p is the enthalpy density, 
e and p are the energy density and pressure, Vi are the components of the fluid velocity, 
and 7 = (1 — |vp)~-^/^ is the Lorentz factor. For spherical bubbles, the only nonvanishing 
component of the stress tensor is T[r) = Trr, and the fluid velocity has only a radial 
component v = Vr- The radial dependence of the enthalpy density w{r) and fluid velocity 
f(r) need to be determined. Gravitational radiation from thin-wall bubbles depends on the 
quantity 
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/ TMr^dr = [ w(r) '"^''^\y dr. (l] 
J J 1 — virr 



The rest of this Section is devoted to evaluating this expression. 

To model a phase transition, we assume that the equation of state of the gas in the 
high-temperature ("symmetric" or "unburnt") phase describes a relativistic gas plus a false- 
vacuum energy density: 

ei= aT^ + e, (2) 
Pi= g^^i^ - (3) 

where e is the false- vacuum energy density (or equivalently, 1/4 of the latent heat). In 
the low-temperature ("broken" or "burnt") phase the equation of state is simply that for a 
relativistic gas: 

62= aT^, (4) 
P2= -aT^. (5) 

Note that Wi = [4:/3)aT^. When a bubble forms in a first-order transition, its interior is 
described by the broken phase equation of state, while its exterior is in the symmetric phase. 
The phase boundary at the bubble wall, the "detonation front", is assumed to be infinitely 
thin. The difference in free energy between the inside and the outside of the bubble creates 
an effective pressure driving the expansion of the bubble. We define the quantity 

a = e/aT^, (6) 

the ratio of vacuum energy to the thermal energy in the symmetric phase; a characterizes 
the strength of the phase transition. The limits a ^ and a ^ oo correspond to very weak 
and very strong first-order phase transitions, respectively. 

In spherical combustion there is no natural length scale, and the hydrodynamic equations 
can be written in terms of ^ = r/t where r is the distance from the center of the bubble and 
t is the time since nucleation. In other words, the velocity and enthalpy-density profiles, 
f(r, t) and io(r, t), are self-similar, being functions of only r/t. The variable ( is then the 
outward velocity of a given point in the bubble profile. As shown by Steinhardt [8], Euler's 
equation and the equations of continuity and conservation of entropy can be combined in the 
case of spherically symmetric flows to yield an equation for the radial velocity as a function 
of^: 



where /x = (^ — 'u)/(l — v^) and 7^ = (1 — v'^)~^ . The enthalpy density satisfles 

1 dw 47V 
w dv 6ci 
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which can be integrated in terms of the velocity profile: 



w{() = Wdexp 



4 r'"d 
3^ 



r 2 , ' 

/ 7 uav 



(9) 



The stress-tensor T[r) can then be obtained from the solutions to these equations with the 
proper boundary conditions. 

Conservation of energy and momentum assure that in the rest frame of the bubble wall, 
/3i, the velocity of fluid in the symmetric phase into the wall, is given by [8,9,12] 



/3i 



{P2 -Pi){e2 +Pi] 
(62 - ei)(ei +P2] 



1/2 



(10) 



and that /32, the velocity of fluid in the broken-symmetry phase away from the wall, is 



/32 



{P2 -Pi){ei +P2] 
(62 - ei)(e2 +pi] 



1/2 



The enthalpy densities on each side of the wall are related by 



l-^l 1-/3: 



(12) 



If Wi [i.e., the temperature outside) and a are given, /3i, /32 and W2 are still undetermined; 
however, once one of the three quantities is given, the other two are determined. 

It has been shown [9] that there are two qualitatively different kinds of combustion. If 
/3i > /32, the transition occurs via "detonation" and the wall propagates at a supersonic 
velocity, i.e., at a velocity larger than c^, the speed of sound; if /3i < /32, the transition 
occurs via "deflagration," and the wall propagates at subsonic velocity. The sound velocity 
is given by dp/de at constant entropy; in general, it is a function of the state variables, but 
in the highly relativistic limit Cg In the remainder of this paper we always take 

this limiting value for the sound velocity. It has also been shown [8] (and will be discussed 
below) that if the transition occurs via detonation, /32 = Cg and so /3i and W2 are given simply 
in terms of a and Wi. However, for deflagrations, /32 is, in general, still undetermined. 

In either case, the fluid velocities (in the rest frame of the wall) in and out of the wall 
are unequal, /3i 7^ /32, so the fluid velocity v must be nonzero somewhere. Moreover, the 
fluid velocity is zero at the center of the bubble (by spherical symmetry) and far away 
from the bubble (in the "rest" frame of the Universe). Thus, deflagration or detonation 
is characterized by a radial fluid velocity proflle, f(r), which satisfles the fluid Eqs. (7) 
and (8) with the appropriate boundary conditions. We now discuss the solution to this 
hydrodynamic problem. 



A. Detonations 

The case of detonations has been discussed in detail by Steinhardt [8], and we review 
the relevant results here. If the transition proceeds via detonation, the unburnt fluid enters 
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the wall at a supersonic velocity. For this reason, there can be no shock preceding the wall, 
so the fluid is at rest outside the bubble wall; i.e., v[() = for ( > where (j, = j3\ is 
the propagation velocity of the wall. Since j3\ > /32, the fluid just behind the detonation 
front is accelerated outward to a velocity (^^ — /32)/(l — (di^2) (this is just the relativistic 
transformation of the velocity from the wall frame to the rest frame of the bubble). As 
shown by Steinhardt, the detonation front is then followed by a rarefaction wave in which 
the velocity proflle v[() falls smoothly to zero at ^ = c^, and remains zero for ( < Cg. 

Steinhardt also showed that detonation solutions to Eq. (7) exist only if /32 = Cg. This is 
the relativistic generalization of the Chapman- Jouget condition for spherical detonations (see 
Ref. [12]). For a general planar detonation [9-11], the value of /32 is not constrained to be Cg. 
Therefore, the detonations in a phase transition in the early Universe, restricted to satisfy 
the Chapman- Jouget condition, are not as general as those considered in some previous 
work [9-11]. We should also point out that the functional form of the velocity and enthalpy- 
density proflles are different in a spherical detonation from those in a planar detonation 
(even with the Chapman- Jouget condition), although they are similar qualitatively. 

Given /32 = Cs, one flnds that the velocity of the detonation front, for a given a is 
simply [8] 

= TT^ • 

In Fig. 1 we plot the velocity of propagation of the detonation front, as a function of a, 
the parameter describing the strength of the transition. The velocity proflle is then given 
by integrating Eq. (7) with the boundary condition ^(^d) = = {(d — Cs)/(1 — (dCg), from 

C = Cd^O C = Cg. 

As ^ ^ (d, dv / d( oo, so Eq. (7) cannot be easily integrated numerically from ( = (d- 
Instead, we write ^ as a function of f, use the relation ([v) ~ + [l/2)("[vd) for v Vd, 
and integrate from some v very close to Vd- Here, 

avd) = -^f^^^'^^- (14) 

is the second derivative of ( with respect to v at the detonation front. The velocity proflles 
for several values of a are displayed in Fig. 2. As shown, v[() is zero for ( < Cg] there is a 
weak discontinuity at ^ = c^, and v increases until ( = (d where dv / d( oo. Also, as a is 
increased, both (d and Vd increase. 

Once the velocity proflle has been determined, the enthalpy-density proflle can be cal- 
culated using Eq. (9). The enthalpy density at the detonation front, Wd = W2, can be deter- 
mined in terms of Wi and a from Eq. (12). Numerically integrating Eq. (9) is straightforward, 
but as the detonation becomes strong, (a ^ 1), t«(^) varies rapidly near the detonation front. 
The quantity /i(^) equals at ^ = c^, increases until some ^ which becomes closer to ^d as 
a is increased, and then rapidly decreases to Cg at the front. One flnds that the region near 
^d where /x is decreasing loosely deflnes a width — which decreases as a is increased — for the 
detonation front, and that t«(^) varies quite rapidly in this region. The enthalpy-density 
proflle, t«(^), divided by lOi, the enthalpy density outside the bubble, is plotted in Fig 3. 
The enthalpy density jumps at the detonation front, then decreases smoothly until i = Cg 
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and maintains a constant value Wq < Wi at the center of the bubble, ( < Cg. For larger a, 
w[() becomes increasingly concentrated near the detonation front. 

In Fig. 4, we plot the stress-energy density T[() = wv'^^f'^. Note that as a ^ 0, all the 
stress-energy becomes concentrated near a thin shell near ^ = c^, while as a is increased, the 
stress-energy becomes dramatically concentrated near the detonation front. The thickness 
of this shell tends to zero in both the strong- and weak-detonation limits and always remains 
negligible compared with the bubble radius; thus a thin-wall approximation to the stress- 
energy tensor of a detonation bubble is valid. 

To gauge the accuracy of our numerical integration, we checked that the energy contained 
in a volume of equivalent size before the bubble was nucleated, -Einitiai = 47r(l + ct)(^/3, was 
equal to the total energy of the bubble: 

^;bubble = ^ tU'l\^+v')dt (15) 

6 Jo 

In all cases, we obtained -Einitiai = -^bubble to within about a percent. 

Although we cannot write analytic expressions for f(^) and analytic expressions 

are easily obtained in the weak-detonation limit, a ^ 0. If a <C 1, then the fluid velocity 
f <C 1 everywhere, and — <C 1. In this case, the small v and ^ — expression of 
Steinhardt [8] can be used to describe the entire velocity proflle: 

^ = c, + ^^(ln^ + l). (16) 



In this limit, Vd = y3a/2, = Cg + y2a/3, and Wd = Wi{l + 2\/2a). As a ^ 0, 7^ ^ 1, 
and /X ^ Cs, so Eq. (9) can be integrated to give 

Wo Wdexp{-4:CsVd) ~ Wd{l - 2V2a). (17) 

The enthalpy densities inside and outside the bubble are equal to lowest order in a, ioq — 
The stress-energy integral, Eq. (1), can be also be integrated analytically. Using d(/dv = 
(2/3) ln(fd/'u), we flnd 

Jcs Jcs 

= wicl / v\dC/dv)dv = {2/27)wiclv^^ 
Jo 

-w-ya^l^. (18) 



2 27 



Since ^ = r/t, the integral over r, Eq. (1), is Eq. (18) times . This should be compared 
with the analogous result for the case of a pure vacuum bubble, Eq. (13) in Ref. [5], which 
in our notation is / T{r)r^dr = Wiat^ /4:. 

In the strong-detonation limit, a ^ 00, both (d and Vd go to unity. Simple analytic 
expressions for w[() and v[() cannot be found in this limit; however, we can flnd a simple 
form for the stress-tensor integral, Eq. (1), using conservation of energy. Equating -Einitiai 
and -Ebubbie, 
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(19) 



In a strong detonation, w[() and 7^ are both strongly peaked at the detonation front, so the 
dominant contribution to the integral comes from values of ( near (d', furthermore, near 
v[() ~ 1, so for a ^ 1, 

Tm'd( ~ w^aej4: ~ a«;i/4, (20) 

which smoothly matches the pure- vacuum result, Eq. (13) in Ref. [5]. 
For arbitrary values of a, we can write 

Z^" T{0(^d( % K{a)wia/4:, (21) 
Jo 

Here, K[a) is an efficiency factor quantifying the fraction of the available vacuum energy, 
or latent heat, that goes into kinetic (rather than thermal) energy of the fluid. Given the 
weak- and strong-detonation limits for the stress-energy integral, Eq. (18) and Eq. (20), and 
the values at some intermediate points that we calculate numerically, we flnd that K[a) can 
be given approximately by 

Aa\—\\ — , (22) 
27 V 2 ' ^ ' 

where A = 0.715. The function /c(q;), along with the numerically calculated value, is plotted 
in Fig. 5. 

No signal precedes a detonation front. Therefore, except for the regions in which the 
bubbles have collided, the dynamics of collision of two (or more) bubbles is simply that of the 
sum of the individual bubbles. This is directly analogous to the case of collision of vacuum 
bubbles, and justifles the use of the envelope approximation for colliding detonation bubbles 
as explained in the following Section. We also mention that the detonation front is stable to 
non-spherical perturbations and therefore remains spherical as it expands [7,12]. It has also 
been recently postulated that although the detonation front is spherical, the fluid behind it 
may undergo a transition to turbulence [13]. We discuss the gravity waves that could result 
from the excitation of a fully developed spectrum of turbulence in the next Section. 
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K{ a 



1 + Aa 



B. Deflagrations 

In Appendix A we present a detailed discussion of the fluid dynamics of spherical de- 
flagration. In contrast to detonations, deflagration fronts propagate at subsonic velocities 
and, as shown in Appendix A, are preceded by a pre-compression shock. However, unlike 
in the detonation case, here T[r) is not concentrated in a thin region (unless the transition 
is weak and (j, happens to be near c^, which we consider unlikely), and the thin- wall ap- 
proximation does not accurately describe the bubble. This makes calculating gravity waves 
from deflagration bubble collisions difficult. However, there are several reasons to believe 
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that the coUision of deflagration bubbles is actually a very weak source of gravity waves. 
Most importantly, the smaller velocities will make deflagrations a much weaker source than 
detonations (the fraction of vacuum energy liberated in gravity waves is proportional to v^). 
In addition, the propagation velocity of the deflagration front is subsonic; therefore, once the 
pre-compression shocks collide, signals can be sent back through the bubble, and there is no 
reason to expect the evolution of two (or more) bubbles to resemble the sum of individual 
bubbles. The spherical shape of the bubble walls is likely to be disrupted shortly after the 
pre-compression shocks collide. Thus, in a deflagration, there will be no large concentration 
of kinetic energy near the bubble walls, so gravity- wave production from the collisions should 
be suppressed. However, we note that deflagration bubbles may be equally as effective as 
detonations at stirring up turbulence, which also leads to generation of gravity waves, as 
discussed in the next Section. 

Actually, the existence of deflagration as a possible mode for a phase transition in the 
early Universe has recently been questioned. It has been argued that a cosmological phase 
transition cannot occur via deflagration because a bubble that begins expanding via defla- 
gration rapidly becomes unstable to detonation due to the existence of hydrodynamic insta- 
bilities [7]. On the other hand, it has also been pointed out that temperature-dependence 
of the propagation velocity of the bubble wall could stabilize a deflagration [11], although it 
seems that this conclusion applies only to very weak transitions. For all of the above reasons, 
and especially the fact that little gravitational radiation is expected from deflagrations, we 
consider only detonations in the following analysis. 

III. GRAVITATIONAL RADIATION 

A. Radiation From Colliding Bubbles 

As in previous work [4,5], we use the linearized-gravity approximation in Minkowski space 
to compute gravity- wave production. In the phase transition considered here, we expect this 
approximation always to be valid; see [4] for a detailed discussion. The fundamental quantity 
for calculating the radiation spectrum is the Fourier transform of the stress-energy tensor: 

r,,(k, Lo) = — dt e^'^* / d^xT,,{-K)e-'''''-^ ; (23) 
zvr Jo J 

we adopt Weinberg's unusual normalization convention [14]. We consistently ignore any pure 
trace pieces of the stress tensor, such as a spatially constant thermal-energy term, as they 
cannot contribute to the production of gravitational radiation. The source here is a number 
of spherical bubbles within a sample volume, each expanding at a given velocity from a 
given nucleation site and time. As a detonation bubble expands, its dynamics until it meets 
another expanding bubble are simple, described by the combustion formalism elaborated in 
the previous Section. Due to its spherical symmetry, a single expanding bubble produces 
no gravity waves. Only after bubble collisions destroy the spherical symmetry of individual 
bubbles is gravitational radiation emitted. In principle, the calculation of gravity waves is 
straightforward: once bubbles are nucleated, simply use the appropriate equations to evolve 
them until the phase transition is complete. For vacuum bubbles, the Klein-Gordon equation 
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is the necessary evolution equation, while thermal bubbles require hydrodynamic equations. 
The stumbling block is the complexity of the bubble configurations once collisions begin. The 
field or fiuid equations in three spatial dimensions require intensive computational resources 
to solve, especially considering the dynamical range in the problem: from the thickness of 
the bubble wall to the Hubble radius. This difficulty prompted the development of the 
envelope approximation [5]. 

In Ref. [4], the full numerical evolution for a pair of vacuum bubbles was performed 
and the resulting gravity-wave emission calculated. The results scale in a simple manner 
with the natural length and energy scales of the problem. In particular, the peak frequency 
of radiation is determined by the size of the bubbles at the end of their evolution, and 
the radiation spectrum varies with the fifth power of this length scale. The results do not 
depend on the smaller-scale structure of the scalar field which develops in the region where 
two bubbles collide. This scaling result suggests that the fine details of the collision region 
are not important to gravity-wave production, but rather that the radiation is dominated 
by the gross features of the evolving bubbles, namely the uncoUided bubble walls. These 
observations prompted the envelope approximation in [5], which consists of treating the 
uncoUided bubble walls as infinitesimally thin energy concentrations and ignoring completely 
the collision regions, in effect considering only the uncoUided "envelope" of the expanding 
bubbles. This approximation turns out to be surprisingly good. In the case of two vacuum 
bubbles, the envelope approximation reproduces the shape and features of the gravity- wave 
spectrum from detailed field evolution, and its amplitude is correct to within about 10%. 
The numerical utility of the approximation is illustrated by a calculation involving nearly 
200 vacuum bubbles nucleated in a sample volume [5], which would be impossible with full 
field evolution even with extensive computational resources. 

As demonstrated in Section II, detonation bubbles satisfy the conditions of the envelope 
approximation. Specifically, the kinetic-energy density is concentrated in a thin shell near 
the bubble wall. In addition, the walls propagate at supersonic velocities, so anything that 
happens in the collision region cannot affect the expansion of the bubble in the uncoUided 
region. On the other hand, defiagrations will not satisfy either condition. First, the energy 
density is not concentrated near the bubble wall; this complicates evaluation of the stress 
tensor, as described below. The most serious problem, however, is that the walls propagate 
at subsonic velocities. This means that the spherical symmetry of the bubble walls can 
be disrupted shortly after the pre-compression shocks collide. Since efficient gravity-wave 
production requires coherent motions of large energy densities, we expect the radiation 
production from colliding defiagration bubbles to be substantially suppressed with respect 
to a detonation of similar strength. 

Using the envelope approximation and ignoring the bubble-collision regions, we can di- 
vide the spatial integration in Eq. (23) into regions, one surrounding each spherical bubble 
centered at the bubble-nucleation site x„. The stress tensor becomes 

r,,(k,a;) = ^ r die'' 
Zir Jo 

where N is the number of bubbles, Sn is the portion of the surface of bubble n that remains 
uncoUided at time t, and the integration variables are chosen independently around each 



N 

E 

.n=l 



JSr, Jo 



dr 



-iwk-x/- 



(24) 
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bubble. If the bubble wall is thin, the exponential can be factored out of the radial integral, 
leaving the r-integral over the profile of the bubble stress tensor independent of the angular 
integral over the uncoUided bubble wall. 

Given the stress-energy tensor, the total energy radiated in gravity waves into a frequency 
interval dco and a solid angle dQ, is [14] 

dE 

= 2G'a;2A,,,,„(k)7;*(k,a;)r,„(k,a;) (25) 
where Aij^im is the projection tensor for gravity waves, 

•'^ij,/m(k) = SilSjm 1\^j\^m8il -\- |kikjk/km ^SijSim -\- |^ijk/km -|- l^/^kikj. (26) 

Contracting with the tensor Aij^im projects out the transverse-traceless piece of the source. 

We model a phase transition by assuming an exponential bubble nucleation rate per unit 
volume [15]: 

r = Toe^\ (27) 

Note that /3 here is unrelated to the velocities /3i and /32 defined in the combustion analysis 
of the previous section. This form is a reasonable ansatz since in general the rate will be 
the exponential of a characteristic nucleation action; keeping the lowest terms in a Taylor 
expansion around the time of the phase transition gives Eq. (27). In general, /3 is expected 
to be of the order 41n(mpi/r)if ~ lOOH for a Hubble rate H [16]. Bubbles are nucleated in 
a sample volume according to this rate. Each bubble expands at a constant velocity until all 
of the sample volume has been converted to the broken phase. The walls of the expanding 
bubbles, treated as thin shells, constitute the stress-energy tensor Tij{pi.,t) in Eq. (24). 

For this form for the nucleation rate, /3~"^ is roughly the duration of the phase transition 
[15], and thus j3~^v is roughly the mean bubble separation (i.e., the bubble size at the end 
of the phase transition). The frequency dependence of the spectrum is set by the time scale 
j3~^ , so the characteristic frequency of the radiation is a; ~ /3. To determine the scaling of 
the amplitude of the radiation spectrum, we note from Eq. (22) that for a single bubble of 
radius R, 



L 



R 1 . . 1 

drr^Tij[r,t) = -R^ K,[a)e^i^j = -R^ K,[a)awi'ki^j (28) 
3 4 



where K[a) is the efficiency factor introduced previously which measures the fraction of 
vacuum energy e converted to bulk motions of the fluid. For vacuum bubbles, k = 1 since 
all of the vacuum energy goes into accelerating the bubble wall. Ignoring for the moment the 
e''^ '^ factors in Eq. (24), Eqs. (25), (28), and (24) imply that for a flxed number of bubbles, 
N, dE/duj oc N{R^Key. (Note that the project ion tensor A contracts with the unit vectors 
in Eq. (28) to form a dimensionless number which depends only on the geometry of the 
problem.) Substituting (3~^v for the length scale gives 

f"^ J oc NG{R^Kef/{Nv^i3-\) oc Gv^K^aw^i3-^ (29) 

duj E^ac 
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where E^^^ ~ NR^e ~ Nev^(3~^ is the total vacuum energy in the sample volume. 

The neglected exponentials correspond to the usual quadrupole approximation, e^^'^ 1. 
Since k-x scales like f, the quadrupole approximation will be valid for small bubble velocities, 
as expected. As v becomes larger, the contribution of the exponentials becomes important, 
and the scaling in Eq. (29) will not hold. In fact, for the case of vacuum bubbles, f = 1, 
the quadrupole approximation overestimates the radiation spectrum by around an order 
of magnitude [4]. Since the quadrupole approximation scales exactly with v^, the actual 
spectrum's amplitude will increase more slowly with v than for larger velocities. Our 
numerical results show that the deviation from scaling begins around f = 0.1; see Fig. 6. 

The radiation spectrum is determined by numerically evaluating the integrals in Eq. (24) 
for the source configuration of many bubbles nucleated in a sample volume. We use tri- 
als with 20 to 30 bubbles because this number is computationally tractable and because 
significantly more bubbles give essentially the same results for the radiation efficiency, as 
demonstrated in Ref. [5]. Thus, for a given value of /3, the physical sample volume is pro- 
portional to v^, insuring that approximately the same number of bubbles will be nucleated 
in the sample volume for any velocity. We have five trial nucleations in a spherical sample 
volume, each with between 17 and 33 bubbles, nucleated randomly according to Eq. (27). 
These are the same nucleation trials used in Ref. [5]. We use the same nucleation trials for 
all bubble-expansion velocities by re-scaling all distances in the v = 1 case by a factor of 
v; using the same nucleation trials minimizes any spectrum differences arising simply from 
geometry of the bubbles. For each trial nucleation and bubble expansion velocity, we cal- 
culate the radiation-energy spectrum in the six directions (ibx, iby, iz), and then average 
over the five trials and six directions to obtain a mean spectrum. These spectra are plotted 
as power per octave for various velocities in Fig. 7. The statistical variation in the mean 
due to the averaging is around 10%. Each spectrum peaks at a characteristic frequency of 
around 2/3 independent of bubble expansion velocity, as expected. In Fig. 6, we plot the 
ratio of energy radiated in gravity waves to the total energy (thermal plus vacuum energy); 
the straight line displays scaling. The departure from scaling as f ^ 1 is clear. The 
solid curve is the analytic fit to the fraction of energy liberated into gravity waves, 

^.O.OT^'flVf^)^^^). (30) 
-Bli \I3) \l+aj \OM + v^) ^ ' 

Note that in the strong-detonation limit, v ^ 1 and a ^ oo, this reduces to the vacuum- 
bubble result of Ref. [3] . 

The radiation spectra in Fig. 7 depend on the parameters f, k, /3, and e = 3ioiq;/4. A 
particular phase transition is characterized by the temperature at which it occurs and its 
latent heat, or equivalently by Wi and a. For detonation bubbles, v and a are related by 
Eq. (13), and k and a by Eq. (22). The parameter /3 describing the bubble-nucleation rate 
will be determined by the effective action for nucleating bubbles. Thus we have assembled 
all the necessary ingredients to calculate the gravity waves produced by a thermal first-order 
phase transition which proceeds via detonation bubbles. 
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B. Radiation From Fully Developed Turbulence 



Injection of energy into the universe will cause turbulence if the Reynolds number of the 
early- Universe plasma is large enough at the time of energy injection. Here we estimate 
the gravity waves produced by a Kolmogoroff spectrum of turbulence, independent of any 
details of the phase transition dynamics. 

The Reynolds number in the early- Universe plasma is very large for length scales L not 
too different than the Hubble radius ~ mpi/T^. Specifically, the Reynolds number 

Re = LV/u ~ 75'^(mpi/r), with L = ^H~^ , the kinematic viscosity v ~ vl, I ~ Ijna ~ 
l/g'^T is the particle mean-free path [g is a typical gauge coupling and T is the plasma 
temperature), and V/v = (bulk flow velocity)/(microscopic velocity) is taken to be of order 
unity. Thus, it is quite reasonable to expect turbulence to develop when the plasma is 
"stirred up" by a phase transition (the critical Reynolds number for the onset of turbulence 
is around 2000), especially if bubble walls are unstable to perturbations and become highly 
nonspherical. 

In the case of fully developed turbulence the distribution of the turbulent kinetic-energy 
density is expected to take the stationary Kolmogoroff form [17], 

j^dpu^ ^ ^_2/3^ ^^^^ 

dk 

which is characterized by a constant flow of turbulent kinetic energy from larger scales to 
smaller scales, 

^ = A%^ = const; (32) 
tl tl dk 

here p is the plasma energy density. The turbulent velocity associated with an eddy of size 
L ~ k~^ , vl, and its lifetime, tl, are related, tl ~ Ljvi,. For the Kolmogoroff spectrum 



vl oc 



L^l^- TL OC L^l\ (33) 



That is, an eddy survives for about a turnover time before it breaks into smaller eddies. (So 
long as the eddy survival time is a scale-independent factor times the eddy turnover time, 
the Kolmogoroff spectrum should develop.) 

On very small scales, k ^ ku, the spectrum is cutoff due to viscous damping of eddies. 
The damping scale ku is the scale on which viscosity diffuses the turbulence as fast as 
the transfer of kinetic energy from larger scales replenishes it: Tdif ~ L'^/i ~ r^; for the 
Kolmogoroff spectrum ku oc On scales k ^ ku, kdpt-arh/ dk oc k~^. 

The Kolmogoroff spectrum is established as turbulence is introduced on some large 
scale — e.g., by the "stirring" of the plasma by expanding bubbles — and is fed down to small 
scales as large eddies break into smaller eddies. It takes of the order of an eddy turnover time 
on the largest length scale to establish the Kolmogoroff spectrum. The stationary spectrum 
of turbulence persists as long as the plasma is being stirred. Once the stirring stops, the 
turbulence dissipates in about a turnover time for the largest length scale. 

Next, let us estimate the amount of gravitational radiation produced by eddies of char- 
acteristic size L. Using the quadrupole formula, Pgw — G[d^Q / dt^Y , and estimating the 
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triple time derivative of the quadrupole moment of a typical eddy as d^Q/dt^ ~ L^pvj^/rL, 
it follows that the volume density of gravitational radiation produced by eddies of size L is 



dpc 
dio 



Lo^f^ ~ Gp'L'Tvl/rl oc Lo-'/' (34) 



where time T is the duration of the turbulence and the characteristic frequency uj t£' 
vl/L ~ vijk. In making this estimate we have made two reasonable assumptions: (i) that 
the quadrupole moment of an eddy varies by order unity on a turnover time; and (ii) that the 
radiation from different eddies adds incoherently. Like the turbulent kinetic energy itself, 
the energy in gravitational radiation achieves its maximum on the largest length scale. 

Finally, let us be more specific. Suppose that the largest length scale on which the 
turbulence is being driven is Lq = /3~"^f, and that the fluid velocities on this length scale are 
Vq [not to be confused with the velocity v of propagation of the bubble wall). Further, we 
assume that the turbulence persists for a time T ^ , corresponding to the length of the 
phase transition. Then we have the following approximate relations: 

VL ^ i^Y'vo; r,^-^ L'I'lTv-'- (35) 



kn ~ [vLollfl^L-^ ~ lb'; ojd ^ T^' ^ vkn. (36) 
It then follows that the spectrum of the energy density in gravity waves is 



LJo ^ Til - /3^"'^o, (38) 

where this spectrum extends from frequency ojq up to ojo. 

Strictly speaking, these expressions are valid only in the regime of nonrelativistic fluid 
velocities, fo <C 1, and likely overestimate the gravity- wave production if applied to a 
stronger transition. For a detonation, the initial fluid velocity Vq can be estimated from 
the fraction of the total energy that goes into kinetic energy of the fluid. Thus, in the 
weak-detonation limit, Vq ~ (/co;)-^/^, and in the strong-detonation limit, Vq ~ 1. For a 
deflagration, the fluid velocity may be estimated by Eq. (Al). 

Our estimate for the gravitational radiation produced in a phase transition should be 
viewed as an absolute, albeit approximate, lower bound. No account was made of the 
radiation emitted by the bubble walls themselves; only that arising from the turbulent 
motion of the plasma that was stirred up by the release of the latent heat was taken into 
account. Further, we wish to emphasize that our analysis and estimates should apply to any 
violent injection of energy on large scales in the early Universe. 
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IV. RELIC GRAVITY WAVES 



To translate the results of the previous section into the potentially observable back- 
ground of gravity waves today, we must propagate the gravity waves forward from the phase 
transition until today. This is simple since the gravity waves are essentially decoupled from 
the rest of the universe. The energy density in gravity waves decreases as and the fre- 

quency of the gravity waves redshifts as , where R is the scale factor. If the universe has 
expanded adiabatically since the phase transition, meaning that the entropy per comoving 
volume S oc R^g[T)T^ remains constant, then the ratio of the scale factor at the transition 
to the scale factor today is given by 

In these expressions, g{T) counts the total number of relativistic degrees of freedom at a 
given temperature, and the star subscript refers to the value of a quantity at the time of the 
phase transition. If we denote the fraction of total energy density in gravity waves at the 
transition as 0.gw* and the characteristic frequency at the transition as then the fraction 
of critical density today 0.gw and characteristic frequency /o today are 



nGw = nGw*[YJ [jfj = 1-67 X lO-^/i-' j Qgw*, (41) 

where h is the current value of the Hubble parameter in units of lOOkmsec"-^ Mpc"-^, and 
we have used the relation 



90m2 



^2 ^ --^ rraa ^ . (43) 



We also define a characteristic amplitude hc{f) produced by stochastic gravity waves around 
frequency / as 

K{f) = 1.3 X 10-''[nGw{f)h'Y^' , (43) 

where 0,Gw{f) is the contribution per frequency octave to the energy density in gravity 
waves [18]. 

Using the results in the previous section, we can describe the gravity waves from bubble 
collisions by 



nG^k-^ . 1.1 X 10-.^ ( f) ( 1 ( 1^^^ (44) 



,2^ r.-e.2iHA f a Vf \ [100 
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« 1.8 X 10-« (— ) (- ] (— ) ( j j ■ («) 

For detonation bubbles, in the weak-transition limit a ^ 0, /c oc -^a, so the amplitude of 
gravity waves is suppressed by a factor of c?!"^ relative to the amplitude in the case of a 
pure-vacuum transition. 

For the case of turbulent mixing, the same analysis applies though our estimates are 
much rougher. We assume that after the phase transition the ratio of the energy density 
in gravitational waves to that in radiation is of the order of ^gw* — {H^^j j3Yva^K^ and 
the spectrum peaks at the frequency 27r/* ~ j3v' Then we have the following 

estimates: 

aowh? ^ 10-= ' vvl {^-^\ , (47) 

/_.2.6xlO-H.„„»-(|-)(^)(^)"', (48) 

M/..>)^5xlO-„5(^)'(i|X)(H^l)'". (49) 

Note that the characteristic amplitude for gravity waves from bubble collisions and from 
turbulence scales in the same way, and our rough estimates indicate that fully-developed 
turbulence is comparable to, and maybe more potent than, bubble collisions in generating 
gravity waves. 

For a particular first-order phase transition, knowledge of the parameters v, j3, k, and a 
suffice to determine the resulting gravity- wave spectrum from bubble collisions. For deto- 
nation bubbles, v and k are functions of a [cf. Figs. (1) and (5). In contrast the time scale 
/3 and the energy scale a are determined entirely by the bubble-nucleation probability. In 
terms of fundamental physical quantities, /3 and a are determined by the effective potential 
for bubble nucleation. Knowledge of the mean bubble separation Lq = j3~^v and the char- 
acteristic fluid velocity Vq suffice to determine the spectrum of gravitational radiation from 
turbulence resulting from the transition. 

As a direct application of our general formalism, we consider the electroweak phase tran- 
sition. This cosmological phase transition has been the focus of much attention recently. If 
the electroweak phase transition was flrst order, then the baryon asymmetry of the Universe 
may have been produced at the electroweak phase transition [19]. Such a transition would 
have produced gravitational radiation; we now use our results to estimate the strength of 
this signal. 

The minimal standard model electroweak phase transition occurs when the SU[2)l x 
U[1)y gauge symmetry is broken to U[1)em- The bubble-nucleation rate and latent heat 
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of the transition follow from the effective potential for the Higgs field (j). In Appendix B, 
we review a general form for the effective potential and its specific realization for a one-loop 
electroweak calculation. We adopt the reference values rrit = 100 GeV for the top mass and 
mn = 60 GeV for the Higgs mass; the end of Appendix B shows how the relevant parameters 
vary with these masses. The transition then occurs at a temperature ^ 104 GeV and 
results in H^/jS = 1.3 X 10-^ a = 1.4 x 10-^ /c = 7.8 x 10-^ and v = c, = 0.57. Then for 
bubble collisions, we get Qh^ ^ 9.8 X 10"^'^ and /ic ~ 1.5 X 10~^^, peaking at a frequency 
around /max ~ 4.1 X 10"'^ Hz. Reasonable changes in the reference values for the Higgs 
and top masses and uncertainties in the accuracy of the one-loop effective potential could 
conceivably change these values by an order of magnitude or more. The weak gravity-wave 
signal that results from the electroweak phase transition is a consequence of the fact that 
the transition in the standard model is very weakly first order, if first order at all. 

Various generalizations of the standard model, particularly enlarged Higgs sectors in 
supersymmetric models, can substantially strengthen the electroweak transition [20]. Other 
more speculative first-order transitions, such as in various GUT theories, may also have 
taken place. We can ask what characteristics must a first-order phase transition possess to 
generate a gravity wave signal which is potentially detectable. For the LIGO facility with 
advanced detectors, the ultimate sensitivity to a stochastic background is an amplitude of 
around 2 X 10~^^ at 100 Hz [18,21]. Requiring the peak frequency of the radiation spectrum 
to fall at 100 Hz, the most sensitive LIGO frequency, gives (/3/i7*)(r*/l GeV) ~ 2 X 10^ by 
Eq. (45). Then for the expected value of /3/i7* ~ 100, Eq. (46) gives ~ 9x10'^^ Ka/{l + a) 
at the peak frequency, making detection by LIGO marginal at best. 

The situation is more promising for a space-based interferometer. Projected capabilities 
of a long baseline interferometer between two satellites are a frequency range from 10~^ to 
10""^ Hz, and a sensitivity down to an amplitude of 10~^^ at 10~^ Hz [18,22]. In this case, 
requiring the peak of the gravity wave spectrum to fall at 10~^ Hz gives (/3 /if*)(T*/l GeV) ~ 
2 X 10'^. Again taking ~ 100, this corresponds to a phase transition temperature of 

20 GeV; the characteristic amplitude of the gravity waves is he — 10~^^Ka/[l + a). This 
background is detectable as long as K,a/[1 + a) ^ 10"'^, a reasonable condition for a strong 
phase transition. These estimates can be made less stringent by noting that the gravity wave 
spectrum for colliding bubbles falls slowly with frequency, and that measuring the gravity 
wave background at a frequency 10 or 100 times higher than the peak frequency only results 
in the amplitude dropping by a factor of a few. We have also not included any gravity waves 
from turbulence, which could give a comparable and independent contribution. A strong 
electroweak phase transition at T = 100 GeV is potentially detectable by a space-based 
interferometer. 

In conclusion, we have calculated the gravitational radiation produced by two potentially 
strong sources during a first-order phase transition: the collision of spherically symmetric 
bubbles, and fully-developed turbulence. Detailed numerical simulation of many colliding 
bubbles leads to a characteristic radiation spectrum which scales with a, k, and /3, param- 
eters related to the latent heat, efficiency, and time scale of the transition respectively; the 
spectrum also depends on the bubble expansion velocity f in a sensible way. Relativistic 
detonation bubbles provide a simple model for bubble dynamics which allows k and v to be 
expressed in terms of a. Likewise, estimates of the radiation spectrum from stationary Kol- 
mogoroff turbulence give similar scalings with these parameters. These estimates indicate 
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that turbulence is likely as potent a source of gravitational radiation as bubble collisions. 
The magnitude of the frequency and amplitude of the resulting gravity-wave stochastic back- 
ground makes detection of a strong phase transition by a future space-based interferometer 
an open possibility, but makes unlikely detection of a first-order phase transition by the 
upcoming LIGO detectors. 
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APPENDIX A: FLUID FLOW IN DEFLAGRATIONS 

Here we present a detailed analysis of deflagration bubbles, analogous to that of deto- 
nations in Sec. II. A. Our aim is to determine the radial- velocity proflle of the deflagration 
bubble. 

We again start with Eq. (7). If we are considering deflagrations, then in the wall frame, 
fluid flows into the discontinuity with a velocity Vi and out of the wall frame with a velocity 
V2 > fi, and both V2,vi < Cg. In the case of spherical deflagration, since the fluid at 
the center of the bubble is at rest, this means that (in the "laboratory" frame) the wall 
propagates at a velocity f2, so the fluid velocity is f = for ^ < i'2- Since V2 > fi, the 
expansion of the gas during combustion exerts a piston effect on the fluid outside the bubble 
and pushes the fluid just outside the bubble with a velocity 

VU = V2) = ^^^^^ = Vo. (Al) 

1 — VIV2 

So in order to determine the radial velocity proflle in a spherical deflagration, we need to 
solve Eq. (7) subject to the boundary condition Eq. (Al). This is straightforward. 

To begin, note that since v,^, (1 — 'u^),7^ > always, dv/d^ < as long as fi < Cg. Since 
fj, < Cs for ( < Cs (the equalities holding only if f = and ( = c^), we know that dv / d( < 
and that v is always decreasing for ( < Cg. The fluid far from the center of the bubble is at 
rest, so for some value of ^ < 1, the fluid velocity v goes to zero. The question is whether 
this occurs for (i) ^ < c^, (ii) ^ = c^, or (iii) ( > Cg. 

If at some value of v ^ 0, then Inv —00, and d(lTiv)/d( —00; however, 
d(lTiv)/d( —00 if and only if the quantity in brackets in the left-hand side of Eq. (7) 
goes to zero (i.e. fi = Cg). Since this does not occur for ( < c^, the fluid velocity v does not 
decrease to zero for ( < Cg. 

Now if we suppose that f ^ at ^ = c^, then we can study Eq. (7) in the limit f <C 1, 
(^ ~ Cs) <C 1, and we flnd that the solution in this case is [8] 

C-Cg='^vln—, v,C-Cg<^l. (A2) 

3 V 
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For V > the right-hand side is always positive, but the left-hand side is negative for ^ < c^, 
so there is no solution to Eq. (7) where the velocity goes to zero at ^ = Cg. 

Therefore, the radial velocity must go to zero for some value of ^ > Cg. Again, iff is to go 
smoothly to zero, then dlnv/d^ — oo, as f ^ 0. It is clear from Eq. (7) that this cannot 
occur for ^ > c^, so a discontinuity must occur, and as we may have guessed for supersonic 
propagation, there must be a shock. Although dlnv/d^ does not diverge as f ^ 0, it does 
go to — oo for some ( > c^; this occurs when fi = Cg [where v = — Cs)/(1 —v()]. So assume 
that this is where the physical discontinuity occurs. Doing so, we find that in the frame 
of the discontinuity, fluid flows into the discontinuity with a velocity /3i = ^ and flows out 
of the discontinuity with a velocity /32 = Cs{^ l/3/3i). In a shock, /3i = 1/3/32 [8], so this 
discontinuity cannot be physical. Therefore, the shock must occur at some value of ( less 
than that at which fi = Cg. 

To flnd the value of ( at which the shock occurs, we again note that in the frame of the 
discontinuity the velocities of the fluid in and out of the discontinuity are /3i = ^ and /32 = /i, 
and then note that in a shock /3i = 1/3/32. This then tells us that the shock occurs when 

Cg Cg 

It is reassuring to note that this occurs for a value of ( smaller than that at which dlnv/d^ 
diverges (determined by fi/cg = 1). 

So, to determine the velocity proflle (and from it the stress-energy tensor) for a spherical 
deflagration bubble, Eq. (7) is integrated subject to the boundary condition, Eq. (Al), until 
fJ-C/cl = 1. At this point there is a shock. As the strength of the transition is increased, Vq 
will increase, and the value of ( at which the shock occurs will increase. This simply means 
that the strength of the pre-compression shock preceding the deflagration front increases as 
the strength of the transition increases. 

Generally, Eq. (7) must be solved numerically, but if the transition is weak, then V2 — Vi 
and fo <C 1. In the limit of small velocities (f <C 1, and as long as ^ — is not too small), 
Eq. (7) becomes 

^ \ ^ _ 2^; 

which can be integrated subject to the boundary condition i'(^o) = '"Oj to give 

According to this solution, near the deflagration front, the radial velocity falls off quadrati- 
cally with radius and then begins to decrease even faster and goes to zero at ^ = c^. Strictly 
speaking, this solution is not valid at ^ = and the radial velocity does not go to zero 
exactly at ^ = c^, but if the transition is indeed weak, the pre-compression shock will be at a 
value of ( just slightly larger than ( = Cg, and Eq. (A5) should provide a good approximation 
to i'(^). In Fig. 8, we plot the fluid velocity as a function of ^ for a rather weak deflagration 
(f2 = 0.1 and Vq = 0.01). We plot the fluid velocity as function of ^ for stronger deflagrations 
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(A3) 



(A4) 



in Fig. 9; the dashed curve illustrates a deflagrations with f2 = 0.1 and Vq = 0.09, and the 
solid curve illustrates the case where V2 = 0.5 and Vq = 0.45. 

The fluid flow in a spherical deflagration is difi'erent from that in a planar deflagration 
[9-11]. In a planar deflagration, the velocity of the fluid between the deflagration front and 
the pre-compression shock is constant. On the other hand, the fluid velocity and enthalpy 
density decrease with increasing ( in spherical deflagration, as we have shown. Therefore, for 
given values of /3i and /32, the pre-compression shock is weaker in a spherical deflagration than 
it would be in a planar deflagration, and in the limit of a weak transition, it is much weaker. 
(Similar conclusions were obtained for non-relativistic deflagrations [23]). Consequently, the 
allowable modes of deflagration in a phase transition in the early Universe may be slightly 
different than those discussed previously [9-11]. 



APPENDIX B: THE EFFECTIVE POTENTIAL FOR BUBBLE NUCLEATION 

Calculation of the gravity waves from a flrst-order phase transition requires two essential 
pieces of information about the transition: the parameters a and /3. These parameters 
characterize the overall properties of the transition and follow from the effective potential 
for bubble nucleation. 



1. A Model Effective Potential 

In a typical flrst-order phase transition, the probability for nucleation of a low- 
temperature phase bubble will be determined by the tunneling action between two vacua of 
an effective potential. To parameterize this effective potential, we consider the general form 

T) = ^-j{T' - T^)f - l«r0^ + Iat^^ (B1) 

where 7, a, and A are arbitrary positive constants and To sets the temperature scale [10,24]. 
This potential possesses two inequivalent minima. The symmetric phase potential minimum 
is always at ^ = where V[(j)) = 0. The broken phase minimum occurs at 




<l> = -iT)^— 1 + (B2) 



where we have deflned 

''-2a' T' ~to-U ^ ^ 

In the second expression for x, we have presumed a quadratic relation between time and 
temperature, tcT^ = I^Tq, valid in a radiation dominated universe at constant entropy. The 
critical temperature Tc at which the free energy of the symmetric and broken phases are 
equal is given by the relation 



21 



At the critical temperature, the energy density of the broken phase first dips below that 
of the symmetric phase; at the temperature To, the symmetric phase becomes unstable. A 
first-order phase transition occurs at a temperature T*, with Tc > T* > Tq. 

To determine the latent heat and vacuum energy associated with the transition, we begin 
with the value of the potential at the broken phase minimum: 



B[T) = -V[v[TlT) 



24A3 



(B6) 



which is the difference in free energy density between the two states of the system. The 
derivative of B is given by 



-7T + 



av[T) 



and latent heat is defined as 



dB 



The vacuum energy associated with the transition is [10] 

e = B{T) - TB'{T). 



(B6) 



(B7) 



(B8) 



To calculate j3 for a given phase transition, the basic quantity we need is V{t) = Ae~^^*\ 
the bubble-nucleation rate per unit volume per unit time. The dimensionful prefactor A is 
expected to be of order but is unimportant for the present calculation. The argument 
in the exponential is the action for nucleating critical bubbles. At high temperatures, this 
action is well-approximated by [25] 



S{t)- 



1237^/2 fT 



a" 



= 13.7aX-^/'x^/'F{x) 
where the function F is defined by 

2.4 



F{x) = 1 + 



1 + 



1 



+ 



0.26 



1 



(B9) 
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This parameterization is accurate to around 1% for < x < 0.95 [25]. 

The nucleation rate is a rapidly increasing function of time near the phase transition, so 
it is sensible to expand the action in a Taylor series about t = t*: [15] 



S{t) ^ S,- /3{t - U] 



'Bir 



/3 



dS_ 
dt 



9A7 1 dS 



I0? tS dx 



> 0. 



:bi2) 



Then the nucleation rate can be rewritten as F = Fq exp /3t as in the previous Section. 
Simple estimates show that the electroweak transition takes place when S ~ 130 [10,25]. 
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2. The Electroweak Case 



The exact parameters of the electroweak symmetry breaking phase transition are not 
yet well known, due both to uncertainties in the standard model {e.g., the top and Higgs 
masses) and to theoretical difficulties in calculating the effective potential, which determines 
the order of the phase transition and the bubble-nucleation rate. For the present calculation, 
we use the one-loop approximation to the finite-temperature effective potential [26] with an 
improved cubic term [25]: where the coefficients are given by 

7 = -^(2m^ + m| + 2m',), (B13a) 



a = — ^(2m^+m|), (B13b) 
ZttVq 



1 
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m 



H 



Svr^ 



-^(2m^ 



(B13c 



= !!i ( 2mt, In ^ + ml In ^ - W In ^] (B13d) 

2v^ IGtt^v^ V "^-^ / 

with vq = 246 GeV, Ina^ ^ 3.51, and In 1.14. 

We adopt the following reference values: W mass mw = 80.6 GeV, Z mass mz = 
91.2 GeV, top mass rrit = 100 GeV, and Higgs mass mn = 60 GeV. With these masses, 
the above coefficients have the values 7 = 0.17, a = 0.019, To = 103.6 GeV, and 
Ao = )^t{T = To) = 0.028. The Higgs self-coupling Ay depends very weakly on T, and 
we will ignore the variation in Ay over the temperature range of interest. 

For the above parameters, x = 0.74 if the phase transition occurs when S = 130. Then 
Eq. (B8) gives e = 0.049r,4 so 

a = 30e/Tr'g^T^ = 1.4 x 10"^ (B14) 

Working out the derivative in Eq. (B12) leads to /3 ~ 400/t*, which gives 

^ = 1.3xl0-^ (B15) 

using the relationship = 0.30mpi/T^gl^' . Since a is so small, Eq. (13) shows that the 
expansion velocity of detonation bubbles is essentially v = Cg = 1 /\/3- Finally, the fraction of 
the vacuum energy which goes into bubble wall kinetic energy is, by Eq. (22), k = 7.8 X 10"'^. 
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FIGURES 



FIG. 1. Velocity of propagation of detonation front as a function of a. 

FIG. 2. Fluid velocity for a detonation as a function oi ^ = r/t for: (a) a = 0.01 (solid curve); 
(b) a = 1.0 (dot-dash curve); and (c) a = 100 (dashed curve). 

FIG. 3. Enthalpy density, ly(^), divided by the enthalpy density wi outside the bubble, for a 
detonation as a function of ^ for: (a) a = 0.01 (solid curve); (b) a = 1.0 (dot-dash curve); and (c) 
a = 100 (dashed curve). 

FIG. 4. Stress-energy density, T'(^) = wv'^j'^, for a detonation as a function of ^ for: (a) 
a = 0.01 (solid curve); (b) a = 1.0 (dot-dash curve); and (c) a = 100 (dashed curve). 

FIG. 5. The fraction k of vacuum energy that goes into kinetic energy of the fluid in a detonation 
as a function of a. The solid line is a numerical calculation; the dashed line is the analytic fit given 
by Eq. (22). 

FIG. 6. The fraction of total energy (within an arbitrary volume) that is radiated into gravity 
waves by colliding bubbles as a function of bubble expansion velocity. 

FIG. 7. The energy per octave radiated in gravity waves for a phase transition with spherical 
bubbles expanding at velocity v, for v = 0.2, v = 0.4, v = 0.6, v = 0.8, and v = 1.0. 

FIG. 8. Fluid velocity for a deflagration as a function oi ^ = r/t for V2 = 0.1 and vq = 0.01. 

FIG. 9. Fluid velocity for a deflagration as a function oi ^ = r/t for V2 = 0.5 and vq = 0.45 
(solid curve) and V2 = 0.1 and vq = 0.09 (dashed curve). 
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